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Abstract 

We develop various complementary concepts and techniques for 
handling quantum fluctuations of Goldstone bosons We emphasise that 
one of the consequences of the masslessness of Goldstone bosons is that 
the longitudinal fluctuations also have a diverging susceptibility char¬ 
acterised by an anomalous dimension {d — 2) in space-time dimensions 
2 < d < 4.In d = 4 these fluctuations diverge logarithmically in the 
infrared region. We show the generality of this phenomenon on the ba¬ 
sis of i). Renormalization group flows, ii). Ward identities, and iii). 
Schwinger-Dyson equations.We also obtain an explicit form for the 
generating functional of one-particle irreducible vertices of the 0(N) 
(non)-hnear cr-models in the leading 1/N approximation.We show 
that this incorporates all infrared behaviour correctly both in linear 
and non-linear a- models. These techniques provide an alternative to 
chiral perturbation theory. Some consequences are discussed briefly. 
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1 Introduction 


The concept of spontaneous breaking of a continuous global symmetry with 
the attendant Goldstone phenomenon of massless scalar excitations has deeply 
influenced many developments in field theory and statistical physics. The ex¬ 
istence of massless modes can already be demonstrated in Landau - Ginzberg 
mean field analysis (or tree level, in the field theory jargon). However, it is 
not sufficient to stop at this level. Precisely because they are massless, the 
quantum/statistical fluctuations of the Goldstone bosons may significantly 
alter the conclusions drawn from mean field theory. Proofs of the Goldstone 
theorem [Q in quantum field theory have traditionally been used only to 
establish the existence of massless poles coupling to the currents of broken 
symmetry. They have not usually focussed on the correlation functions. Most 
of the applications in quantum field theory have been based on the under¬ 
standing at tree level. Notable exceptions are the soft pion theorems P| 
combining PGAG with current algebra.These are mostly applied to extract 
the leading infrared behaviour. A general analysis of these issues based es¬ 
sentially only on symmetry arguments was given some time ago by Weinberg 
[|]. In that article he also proposed methods to go beyond the leading order 
results. He also applied renormalization group methods to understand some 
universal features of these results. Quite recently, some of Weinberg’s ideas 
have been extensively used in the context of the so called “chiral perturbation 
theory” §. 

In contrast, in statistical physics, powerful techniques have been devel¬ 
oped for handling the fluctuations in this context. These are of importance 
for quantum field theory also.Some examples are: i)the finite temperature 
chiral symmetry restoring transition in QGD, where Goldstone phenomenon 
in three Eucilidean dimensions is of relevance |^; ii)QGD at temperatures 
below the chiral transition temperature as well as zero temperature QGD, 
where Goldstone phenomenon in 4 euclidean dimensions is of relevance,and 
iii)the problems that are handled by chiral perturbation theory at present, 
i.e. effects of pion loops on various processes 0. 

In this paper we point out various complimentary concepts and techniques 
for handling the effects due to quantum fluctuations of Goldstone bosons. 
We develop techniques for computing correlation functions away from the 
infrared regime.We show that there are certain universal features, valid even 
in four space-time dimensions.Our approach provides an alternate procedure 
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for chiral perturbation theory which is closer in its spirit to what has been 
proposed by Weinberg on the basis of RG equations. This aspect is elaborated 
upon in section 7. 

2 Infrared dynamics of pions and sigma 

Often the paradigm in quantum held theory is what is seen from a tree level 
analysis of the 0(N) linear and non-linear cr-models. The linear cr-model is 
given by the action [Q, 

s = / - cr) (1) 

$ = transforms as the fundamental representation of 

0(N). If > 0, the minimum of the classical potential is at = C^. 
Therefore the VEV < <h > is no longer zero in a tree level analysis.Choosing 
the vacuum < >= G, < >= 0, i = 1, 2, ...(A^—1), among the equivalent 

vacua, we get, 


S = J + ^{d^a{x)f -UC^a'^ix) 

-UCa{x){a‘^{x) + vf^(x)) - + 7f^(a;))^ (2) 

where a = ^n— < >,7? = {$^,7 = l,2,...(iV — 1)}. Therefore at tree 

level, we have {N — 1) massless pions and a massive a with a mass v2UC^. 

Precisely because the Goldstone bosons are massless, loop corrections 
may drastically alter the above picture. Now it is well known that infrared 
divergences due to massless huctuations drastically alter the conclusions of 
a semi-classical analysis, especially if the spacetime dimension d is less than 
four P .Already, way back in 1940, Holstein and Primakoff 0 argued that 
not only the transverse mode (i.e. tt) but also the longitudinal mode (cr) is 
soft in a quantum ferromagnet in three dimensions. Their arguments can be 
understood,as follows, in terms of the non-linear a -model. 

The non-linear a- model is relevant to the infrared dynamics of Goldstone 
bosons [§] .A simple heuristic argument is that if a is massive, then for 
the study of the infrared dynamics of the pions, we may use an effective 
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Lagrangian obtained by integrating over the heavy a degree of freedom. Also 
for leading infrared behaviour we may retain only the lowest (two) derivative 
terms in the effective action.The result is simply the Heisenberg model, 

= ^ / (3) 

with the constraint, (j)‘^{x) = 1. Statistical physics in d space dimensions is 
simply the Euclidean quantum held theory in d space-time dimensions. In 
statistical physics t is proportional to the temperature, whereas in quantum 
held theory it is related to the pion decay constant, t = 2/“^(in what fol¬ 
lows, we set /tt = 1). If the spins are ordered in the direction of (f)N, then 
eliminating this held using the constraint, the action can be written entirely 
in terms of if = {0j, i = 1,2, ...{N — 1)}. They describe the spin waves or the 
transverse huctuations.At the tree level for the vr propagator we get. 


K{k) = k 


-2 


( 4 ) 


We may also consider the correlations for longitudinal huctuations i.e. 

of the held a = cpN I — vf^ : —1(: refers to any ordering prescription 

that makes products of helds at the same point meaningful). This means, 

a : -|-. We may expect the leading infrared behaviour to 

come from the hrst term. This would mean that so far as the leading infrared 

,2 


behaviour is concerned. 
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a 


TT^ ;. If we ignore self-interactions of 


vr, and use the above propagator for vr we get,for the a propagator. 


A,(a:) = Al{x) 


X 


\A-2d 


This gives, atleast for d < 4,in momentum space. 


A.(fc) =1 k 


( 5 ) 

( 6 ) 


Ignoring the self interactions of tt’s will later be seen to be justihed because 
these are soft in the infrared region- a consequence of the soft pion theo- 


k 


1-1 


to be contrasted 


rems.The result is that in d = 3 the cr-propagator is 
with the vr-propagator.This means that a also has long range correlations. 
This appears to negate the naive deduction of the non-linear a model ob¬ 
tained by integrating over the heavy a held in the linear a model |T^ . 
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One could dismiss these results as artifacts of a 1-loop calculation. In fact 
the infrared divergences grow with the loop order. Therefore it is impera¬ 
tive to tackle this issue in a framework that transcends loop expansion.There 
is an exactly solvable model, the Berlin-Kac model which goes beyond 
loop-expansion and again gives precisely the same diverging longitudinal sus¬ 
ceptibility and also a massless a in d < 4. Again it is not clear whether this 
is an artifact of the model. 

In order to get a handle on the infrared dynamics of the a held, Pokrovsky 
and Patashinky proposed a ’’Conserved Modulus Principle” (CMP).This 
principle holds that the huctuations are such as to maintain the modulus of 
the N- component order parameter. In mathematical terms CMP states that 


2 < (j) > 6(j)L + 6(j)^ = 0 (7) 

Now the expectation value of is given by the Ornstein- Zernike form: 


= J 




q'^ + H/ <(j)> 


( 8 ) 


It is convenient to use susceptibility to describe the infrared behaviour.The 
longitudinal and the transverse susceptibility in presence of an external mag¬ 
netic held H = (0, 0,..., 0, H) are given by , 


XL = limk^oA^{k^), xt = limk^oA^{k^) 


Using Eqns. and ||, we get, for d < 4, 

(^)d/2-2 


d{6(j)L) 

XL = - ~ 


dH 


Xt 


H 


-1 


(9) 


( 10 ) 


Thus the longitudinal suceptibility also diverges. Patashinsky and Pokrovski 
have ohered a proof for CMP in ||^|. They express the change in the ther¬ 
modynamic potential U computed to quartic order in the huctuations as 


5U = l/2[ 




( 11 ) 


Xt 

and interpret this as implying CMP. Though to quartic order the correct 
form of this variation, which dihers from the one given above (and which can 
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be verified using the Ward identities discussed later), is 


5U 


1 

2xt 




, 1 Mr 

402 ^ 2xl 402 


( 12 ) 


the fact that CMP is equivalent to the correct statement that the extreme 
IR behaviour is governed by the non-linear cr-model indicates the possibility 
of a better proof of CMP. 


The new era for understanding the role of fluctuations of Goldstone 
bosons started with the renormalization group techniques. Brezin, Wallace 
and Wilson computed the equation of state for 0(N) a-model using the 
e-expansion.They found important differences from the corresponding calcu¬ 
lation for the N=1 case (i.e. the Ising model).This could be traced to the 
infrared singularities due to the Goldstone bosons that are present every¬ 
where in the coexistence region for N > 1. In particular it meant that the 
conclusions of the tree level analysis could not be right (atleast for d < 4). 

At the critical point {T = Tc, H = 0), a and tt are on the same footing. 
Both their susceptibilities diverge with vanishing magnetic field as, 


Xt = Xl^ Hs 1 (13) 

governed by the critical exponent 6. The question of interest is the behaviour 
for T < Tc as H ^ O.The assumption that the singularities due to Goldstone 
bosons exponentiate in the e-expansion suggested that. 


Xt 


H 


-1 


XL 




(14) 


Holstein -Primakoff |^, Berlin-Kac [^, and Patashinski-Pokrovski |]^ 
models have all made the same prediction. Such a result has also been proved 
rigorously for 0(2) case using correlation inequalities in Ref. . 

A deeper understanding can be obtained by considering the renormaliza¬ 
tion group flows shown in the figure below: At the critical point, T = Tc, H = 
0, the infrared dynamics is governed by the 0(N) non-Gaussian fixed point 
MQn when 2 < d < 4.Everywhere on the critical surface the renormaliza¬ 
tion group (RG) trajectories flow from the the 0(N) Gaussian fixed point 
(which is infrared unstable) to the 0(N) non-Gaussian fixed point MQn- For 
T > Tc we have the symmetric phase with massive (and degenerate) a and 
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Figure 1 

Figure 1: The fixed point structure of the 0(N) sigma model; the arrows 
indicate infrared flows 


vf.The infrared behaviour is governed by the trivial fixed point T which cor¬ 
responds to zero correlation length.T < T^. presents a new feature compared 
to N=1 case (i.e.the Ising model). The theory has long range correlations for 
all T < Tc, if = 0 as a consequence of the Goldstone phenomenon.Hence the 
infrared behaviour must be determined by a fixed point of infinite correlation 
length.What is this fixed point? 

T — Tc, the deviation from the critical temperature in statistical mechan¬ 
ics, is to be identified with the square of the renormalized mass IQ] in the 
linear a-model. In particular they are parameters with positive mass dimen¬ 
sionality. Under RG flow they necessarily scale.Thus, once T — Tc < 0, the 
flow is towards T = O.This means that we are interested in the zero temper¬ 
ature fixed point. An indication as to the nature of this fixed point came 
from a study of 2 -1- e-expansion of 0(N) non-linear cx-model by Brezin and 
Zinn-Justin jl^ . 

The 0(N) non-linear a- model in d=2 is renormalizable. Even though 
the renormalized perturbation is in terms of the spin wave modes, Mermin- 
Wagner theorem [Q requires that the system is never in a phase with sponta¬ 
neously broken 0(N) symmetry. A hint of this is contained in the renormal¬ 
ized perturbation theory: the theory is asymptotically free.The coupling con¬ 
stant grows with the distance so that the long distance structure of the theory 
could be drastically different from the perturbation theory spectrum.Indeed 
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the spectrum is believed be that of massive and degenerate a and tt’s in the 
fundamental representation of 0(N). 

A consequence of the asymptotic freedom in d = 2 is that in d = 2 + e, 
the /3-function /3(t) has a non-trivial zero (say at t = tc) in addition to 
the trivial zero, f = 0, corresponding to the theory of (N-1) free scalars.At 
t = tc the theory is scale invariant and corresponds to the critical theory 
with N massless scalars and with anomalous dimensions ( denoted by MQ*^ 
).lt should be noted in Fig. 1 that the hxed point NQn of the linear 0{N)- 
model is identihed with the hxed point AfQ% of the non-linear model. This 
follows essentially from a remark made in p[ that at MQn the correlation 
functions of the linear model are identical to those of the non-linear model 
where in the latter one includes correlation functions of composite operators 
also. The coincidence of these two hxed points can also be established in 
the 1/iV-expansion scheme |^. For t > tc, the coupling constant grows 


with the distance, and as in d = 2, the theory is in the unbroken phase, 
f < tc is of interest to us.In this case the infrared how is towards t = 0. The 
infrared behaviour is of O(N-l) Gaussian theory (This hxed point is denoted 
by Qn-i in Fig.l.) and the ultraviolet behaviour is of the critical theory with 
anomalous dimensions. 

The fact that the infrared behaviour for t < tc is governed by the O(N-l) 
Gaussian hxed point gives exhaustive information on the ehects of quan¬ 
tum huctuations of the Goldstone bosons.lt means that the pions decouple 
from each other at low momenta, and their propagators behave as the free 
propagators for small momenta.lt also allows us to obtain the non-leading 
IR-behaviour as well as the ehects of an external magnetic held.In particular 
we now understand why various heuristic considerations all gave the correct 
answer: it was the fortuitious behaviour of a free massless theory.Anything 
more complicated would have been much more difficult to handle. 

Upto what dimension d is the analysis of d = 2 -|- e calculation correct?ln 
particular is it applicable to d=4?The case d=4 is somewhat diherent be¬ 
cause of the absence of the non-Gaussian hxed point A/'^Ar.At the critical 
temperature the infrared behaviour is believed to be given by free massless 
theory of N scalars as the beta function is positive for small couplings; there 
is no evidence for a non-Gaussian hxed point.This appears to exclude the 
validity of the 2 -|- e analysis to this case. We need an alternate technique to 
handle this situation.This is provided ,for example,by the 1/N expansion. 

The 1/N expansion provides results which agree with the above analysis 
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for all d.(The Berlin-Kac model is a variant of this technique.) It yields 
the non-Gaussian hxed point for 2 < d < 4 and shows its absence for d = 4.In 
either case it shows that the leading infrared behaviour anywhere in the 
broken phase is given by N-1 free pions.The only handicap of this technique is 
that it does not provide non-anomalous dimension for the non-Gaussian hxed 
point MQn in the leading N approximation.However even this is overcome 
once one takes higher order (in 1/N) corrections. There are also indications 
that the technique gives qualitatively correct and quantitatively reasonable 
results for all N > 2. 

The fact that even in d=4 the pion dynamics is governed by the 0(N- 
1) Gaussian hxed point Qn-i in the infrared has important implications.In 
particular the propagator of a (to be interpreted as the 0(N) partner of the 
pious) has an universal logarithmic singularity in the infrared region, 


( 16 ) 

where p is the renormalisation scale. This follows from the same arguments 
as for d < 4, in particular because a vr^ : in the infrared.This feature is 
present whether cx-propagator has a pole at some 7 ^ 0 or not. The precise 
location of the pole is a non-universal feature depending on the microscopic 
dynamics. It could have interesting implications for nuclear forces at low 
energies in the ’sigma’-channel (i.e., spin J=0,isospin 1=0) as will be discussed 
elsewhere. 

In Sec. 3 and 4 we develop a new formalism for N ^ 00 which will allow 
for explicit calculation of all correlation functions even away from the infrared 
regime.We will explicitly demonstrate all the properties mentioned above.Our 
formalism gives the same results in the leading infrared limit whether one 
starts with the linear or the non-linear a- model. 

That the pions behave like free massless bosons in the infrared limit, in 
any d > 2 , is such a general phenomenon that there should be an easier way 
of understanding it. A careful analysis of the proof of the Goldstone theorem 
along with current algebra and PGAG provide one such non-perturbative 
argument.Another is the treatment based on Ward identities given in sec 
5.Low energy theorems based on PGAG and current algebra do imply that 
pions decouple from each other and with other matter at low momenta, but 
they do not directly imply that the pion has a canonical dimension in the 
infrared. Gurrent algebra requires that the current coupling to the pion has 





canonical dimension d — 1. However, the canonical dimension (for the k ^ 0 
behavionr of the pion) is (d — 2)/2. This only means that the PCAC relating 
the two needs a dimensionfnl parameter, bnt does not directly relate the two 
dimensions. 

Consider the time ordered correlation fnnction of a cnrrent in a broken 
direction and the corresponding vr held: 

A”W=<0|T(j‘(i),;r‘(0))|0> (16) 

We get, 

J d^xd^^A^^ix) =< 0 I [Q\7r*(0)] | 0 0 (17) 

since [Q*, vr*(0)] = ct( 0) has a non-zero expectation valne. This means for the 
Fonrier transform of this correlation fnnction, 

hmg^oq^A;{q) ^ 0 (18) 

This reqnires that in the spectral representation of the correlation fnnc¬ 
tion there is necessarily a d(m^) contribntion implying a simple pole at zero 
mass.Any other singnlarity, snch as a branch point will not satisfy the above 
eqnation. 

The fact that the a channel also has an infrared singnlarity is eqnally 
general and shonld be nnderstandable in a general way. In Sec. 5 and 6 we 
provide two proofs which rely essentially on symmetry argnments. The hrst 
nses Ward identities, and the second, the Schwinger- Dyson eqnations. 

3 r in ^ (X) limit 

It was argued in Sec. 2 that the leading term in 1/A^ expansion retains all im¬ 
portant effects of the Goldstone boson fluctuations. Analytic calculations are 
possible in this limit and therefore we can analyse such effects in detail.There 
is an extensive literature on this technique |§[.In this section we provide a 
new arsenal to this technique.We dehne a generalized generating functional 
of one-particle irreducible (IPI) vertices and obtain an explicit expression for 
it in the large N limit. This is made possible by involving an auxiliary held. 
Our method provides the easiest way of extracting all correlation functions 
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without doing a separate calculation for each |^. We present results for 
both linear and non-linear cr-models. 

W[H]^ the generating functional of connected Green’s functions, is given 


by, 





(19) 

where. 




H.^ = [ d^xH{x)M{x) 

(20) 


Connected Green’s functions of are obtained from various functional deriva¬ 
tives of r[$], the generating functional of IPI vertices, is obtained 

from by a Legendre transformation. 


T[$] = W[H]-H.^ (21) 

where. 




6W[H] 

6Hi{x) 


( 22 ) 


The vacuum expectation value (VEV) of <hj is obtained as the value at 
which r is minimized.The second functional derivative about this minimum 
gives the inverse propagator. Further, various functional derivatives of T 
about this minimum give the IPI vertices. The connected Green’s functions 
can be extracted by building only tree diagrams using these propagators and 
vertices. 

In the linear a-model, introduce an auxiliary held A to make the exponent 
bilinear in <h: 


^W[H] _ j X>XD^e^ d‘^x(^(d^S(x))^-^A(x)($^(x)-C^)+^\^(x)+H(x)-S(x) ^23^ 

Now an integration over each <hj,i = l,2....,iV, gives a factor. 

Therefore, 






(24) 
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where, 


S/{x,y) = {-d'^ - X{x))5{x-y) (25) 

On defining, 

NU = u,C^ = Nc^, Hi = VNhi, W[H] = Nw[h] (26) 

We get, 

^iNw[h] _ J rf‘* 2 :(^A 2 + ic 2 A))+ 5 TWnV- 5 /i.V-l.h] ^ 27 ) 

We are interested in the limit — >■ cx), f/ —> 0 holding NU = u fixed. 
With the variables redefined as in Eqn. ^ it is clear that this limit is 
dominated by the saddle point. The fiuctuations about the saddle point give 
contributions which are non-leading in N. Thus the leading contribution is, 

Woc[h] = [ d^x{^X\x) + ^c 2 A(a;)) + ^Tr/nV - (28) 

t/ tJj ^ ^ ^ 

where A(a;), the saddle point, is a functional of h{x) as given by, 

-A(a:) + - iV'Ux, x) - (h.V-^){x).{V~\h)(x) = 0 (29) 

u 

We may obtain r[<h] from taoo[^] by a Legendre transform as described 
earlier. But this form is not easily amenable to further analysis. Instead we 
note that the saddle point equation may be rewritten as 

—X{x) + c‘^ — i'V~^{x,x)=(l)‘^{x) (30) 

where, 

X7(j)i{x) = -hi{x) (31) 

Notice that these equations can be obtained from the following functional of 
0 and A : 


roo[0, A] = J d'^xi^X^ix) + ^c^A(a;) + 

— ^X{x)^{x)) + ^TrlnV 


(32) 
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by taking functional derivatives: 


- hi{x) 


<5roo[0,A] 

5(t)i{x) 


—k(x) 


6T^[$,X] 

6X{x) 


(33) 


Here we have temporarily added a ’source’ k for A.When restricted to corre¬ 
lation functions with only (f)i on external legs, k can be set to zero, yielding 
eqn 30. These equations imply that Woo[h, n] and Tooi^, A] related by 
a functional Legendre transformation in both arguments. This means that 
the connected Green’s functions of 0 and in fact also of 0^ (i.e. functional 
derivatives of Woo[h,K\ wrt both h and k) can be obtained from tree dia¬ 
grams built from the ’action’ roo[0, A].This connection is simply an algebraic 
consequence of the Legendre transformation. We have thus the option of 
constructing Green’s functions involving the composite operator : 0^ : also 
by having A as external legs. 

By introducing an auxiliary held A we have obtained an explicit form of 
the generating functional of IPI vertices which is far more amenable to anal¬ 
ysis than if we had obtained a r(0) by substituting the saddle point solution 
for A .Moreover we hnd that the properties of the phase with spontaneously 
broken symmetry can be understood in a simple way in terms of this aux¬ 
iliary held. Ofcourse, now A also appears in propagators and vertices. The 
usual generating function T [0] can however be recovered by using k = 0 from 
Eq.^ to replace A as a functional of 0 in roo[0, A]. 

There is hardly any change in our calculations (and results) when passing 
over to the non-linear cr-model. Now we have. 


Using the Fourier representation for the functional 6 function. 


In this form the only diherence with the linear cr-model is that a quadratic 
term in A is missing in the exponent.Equivalently, the non-linear cr-model is 
simply the special case U —> cx) of the linear cr-model. Thus all results for the 
non-linear cr-model can be extracted by taking this limit in our formulae.We 
hnd that it is straightforward to take this limit and the results are not very 
diherent. 
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We have been cavalier in our manipulations with the functional inte¬ 
grals.For instance, the Trln term has to be carefully dehned.We may choose 
a Pauli-Villars or a lattice cut-off for dehning this term.If we consider a func¬ 
tional Taylor expansion of the Tr In term about a constant value of A, for 
2 < d < 4, only the term linear in A is sensitive to the cut-off.We may there¬ 
fore obtain a renormalized roo[0, A] by making a specihc cut-off independent 
choice for this term. In d = 4, in addition, the term quadratic in A is also 
divergent and requires a renormalization.For the sake of uniformity, we shall 
introduce local counterterms aX + hX^ for both d = 3 and d = 4. In the 
d = 3 case, this amounts to a hnite renormalisation of u. Because of the 
presence of massless modes, renormalisation should be so performed that it 
does not introduce spurious infra-red divergences. Hence we choose to de- 
hne renormalised quantities at A = /i(say). Then the renormalised effective 


action is 



1 

2 


X{x)^‘^{x)) + ^Tr'^lnW 


(36) 


Tr'^ means that the local part of the hrst two terms in the functional expan¬ 
sion around A = /i are to be ignored.Now u and c are renormalised quantities. 
Thus the model is renormalizable in the usual sense.lt is also seen that the 
wave function renormalization for A, equivalently the composite operator 0^, 
is hnite in the iV —> cx) limit.(For a discussion of renormalizability of the 1/N 
series see Ref. [^). 

The discussion of renormalizability is equally valid for the the non-linear 
(T-model. The only difference is that the bare A^(x) term is not present in the 
action.The quantum effects introduce such a term which is hnite in 2 < d < 4 
and diverges logarithmically with the cut-oh in d = 4.Thus the non-linear 
(T-model is not renormalizable in this sense in d = 4 even in the iV —>■ cx) 
limit. 

4 Correlation functions in N ^ oo limit 

In this section, we extract the infrared behaviour of the correlation functions, 
in particular, in the phase with SBS. 
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First we need to obtain the VEV’s Vi and Aq of the helds (j)i and A.For 
uniform (in particular zero) helds h (and k), the VEV’s are space- time 
independent.Therefore we get 


XqV = h 

d'^k , 1 




Aq — /i 


(37) 

(27r)‘^— Ao k'^ — fi {k"^ — ^ ^ 

The analysis below closely follows Ref. [^. We hrst consider the situation 
in the absence of the external magnetic held, h = O.Now AqU = 0 and so there 
are three cases: 

i) Aq = 0 and n = 0 corresponds to the critical point.This is reached for 
the critical value c = Ccr where, 


+ f— _i_ 

"" M ^ 7 (27r)'^F(F-/i)2 


ii) Ao 7 ^ 0 and n = 0 is the symmetric phase. Aq is given by. 


(39) 


An(- i 


d^k -|- AqA:^ — 2pfc^ 


r) = (d - c^) 


u J {2 tiY k‘^{kP‘— XojikP'— jJiY' 
where we have used Eqn.|3^ 

iii) Ao = 0, n 7 ^ 0 is the phase with SBS.The VEV v is given by. 


(40) 


= (c^ - d) (41) 

as seen by using Eqn.|39| 

For any Y > we have a solution for a real v and this corresponds to a 

phase with SBS. On the other hand for any Y < we have a solution for 

Ao in Eqn.|40|, corresponding to the symmetric phase. 

In terms of the huctuations 

A'(a:) = A(a:) — Aq (42) 
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we have for the symmetric phase, 

roo[0, A] = J - ^\o(j)‘^{x) - ^\'{x)(j)‘^{x)) 

+^J J d‘^xd‘^yX'{x)I{x-y,\o)\'{y) 

+lTr'^H-d^ - Ao - A') (43) 

Where I{x — y, Aq) is the fourier transform of 

^(^'’Ao) =*/'^''^[(A;2_^)2 “ _ Xo){{k + p)^ - Xo)^ ^ U 

All terms linear in the fluctuations drop out because we are expanding 
about the extremum value.A q is seen to be the {massY of the 0 field. 

For the phase with SBS we choose utv = 7 ^ 0, Uj = 0, i = 1, ...(iV — 1). 

In terms of the fluctuations, 

a{x) = (j)N{x) — V (45) 

we get, 

roo[<h,a. A] = J d‘^x{^{d^(j){x)f-vX{x)a{x)-^X{x){a‘^{x)+Tf‘^{x)) 

0)A(l/) + i]^Tr'^ln{-d‘^ - X) 

(46) 

TT is massless as is to be expected. The mixing of a with the auxilary held A 
is a crucial feature of this phase. This mixing governs the infrared behaviour 
of the correlation functions as seen below. 

The terms quadratic in a and A may be written in momentum space as, 

i{-X{—k)I{k‘^, 0)A(A;) — vX{—k)a{k) + -a{—k)k‘^a{k)) (47) 

Now, 

li^) = J(A;^ 0) = ^ - 2<d<4, 

Zi (X 

~ — +/n(/i//c^), d = 4 (48) 
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We get the propagators by inverting the matrix of the quadratic form in 

Eqn.^. 


< a{—k)a{k) > 

< a{—k)X{k) > 

< X{-k)X{k) > 


P - r;7/(F) 
iv 

P/(F) - 
ik^ 




(49) 


We are now in a position to extract the infrared behavior of various 
correlation functions. The propagators have the following infrared behaviour: 


< 7ii{—k)nj{k) > 


< a{—k)a{k) > 

< a{—k)X{k) > 

< X{-k)X{k) > 


rsj 


rsj 


rsj 


rsj 


•> 


I /c 1^ 
i 

n’ 

ik'^ 



(60) 


Thus the leading N approximation reproduces the correct infrared be¬ 
haviour for the 71 and the a propagators, viz., tt has canonical k~^ behaviour, 
whereas a has an anomalous dimension (d — 2), precisely as if a vr^ :.We 
now demonstrate explicitly that the other correlation functions also have the 
infrared behaviour required by the soft pion theorems. 

The interactions of tt and a are only mediated through the auxiliary 
held A in this formalism, and the A propagator has an unusual k"^ infrared 
behaviour.This is the essential reason why tt’s and a ’s decouple from each 
other in the infrared, and the results of soft pion theorems are valid. 

For the leading IR-behaviour, the multi-A vertices are irrelevant. The 
effective action is then quadratic in A and consequently, the A-held can be 
exactly integrated out.This leads to the following effective action to describe 
the leading IR behaviour: 


YIR 


d-‘k[^\mr-+ 


\{}^ + 2va){k)\^ 

4/(F) 


(51) 
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This summarises all the essential results of this paper. 

Thus the large N approximation provides a valuable non- perturbative 
tool for obtaining effects due to quantum fluctuations of Goldstone bosons.lt 
incorporates the infrared effects correctly.lt also incorporates the ultraviolet 
behavior correctly: for 2 < d < 4, the ultraviolet behaviour is of the non- 
Gaussian fixed point, and for d=4, it is of the 0(N) Gaussian fixed point. 
It has an explicit a field, so that effects due to its infrared behaviour on 
other matter can be investigated. Also the technique provides a way of going 
beyond the leading infrared region.lt provides a detailed structure for the 
propagators of both n and a with correct asymptotic behaviours. 

We now consider the situation with a uniform external magnetic field 

h = (0, 0,. ,h). Now, AqU = h implies that neither Aq nor v is zero.This is 

the choice that removes terms linear in the shifted field a one gets essentially 
the same Too as the one we got for the symmetric phase bnt with a crucial 
additional term —vX'{x)a{x) and with the proviso that elsewhere 0^ now 
stands for + vf^. 

Thus TT has a mass V^o- has a mixing with A as in the SBS case.The 
aX propagators may be computed as before. 

The equation governing Aq is now. 


-Ao = {v^ - Sc^) + iAo 
u 


dX^k , 1 

(27r)'^ — Aq) 


1 


{k"^ — /i)2 


(52) 


with — c^j.. By an analysis of these eqnations one obtains a mass 

that scales as for the pions and as for the sigma. 

For applications to statistical mechanics, one should use the expressions 
derived here continned to the relevant nnmber of enclidean dimensions. 


5 Use of Ward identities 

In this section we demonstrate that the xl diverges for 2 < d < 4 ,that a ~ 
in the infrared regime and the softness of mnltipion scattering amplitndes 
using Ward identities. 

Ghiral Ward identities merely reflect the nnderlying symmetry strncture 
of the theory and are usefnl both in the symmetric phase where they relate 
correlation fnnctions of the same type ( n-point fnnctions with fixed n) as 
well as in the ordered phase where they relate different n-point fnnctions in a 
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non-trivial way. Note that in the generating function both the measure V(j), 
and the action S'(0), are invariant under 0{N) rotations generated by T^dhe 
source term is also invariant under a simultaneous rotation of both the source 
H and the 0 held.Considering an inhnitesimal transformation, this implies 
for r, 


when an external held H is present.This means that T has only 0(N) in¬ 
variant combinations such as 0^, (0^0)^ etc., except for one additional term 
/ depending on the external held. 

By successive diherentiations of this identity w.r.t. the helds a, vf, one 
can generate various forms of Ward Identities.In particular, an identity that 
will prove most useful later on is 

r<g(p) - rg(p) =<a> r.,„(p,o, -p) (54) 

where < cr > is the expectation value of the order parameter. 

To explore the infrared properties, we consider the thermodynamic poten¬ 
tial, which is called the the exact ehective potential in quantum held theory. 
We include an external magnetic held to serve as an infrared regulator and 
consider the singular behaviour as this held goes to zero. The various zero 
momentum n-point functions are obtained by functionally diherentiating this 
potential w.r.t. constant helds and evaluating the derivatives at the minima 
of the potential. In the ordered phase this minima occurs at a non-vanishing 
value < (T >= n as the magnetic held is reduced to zero. 

As a consequence of the Ward identity, Eqn.^ the ehective potential has 


the form. 

17(0) = /(0-0)+#-0 

(55) 

Since 

U, = ^ = 20,/' - H, 

0(pi 

(56) 

the VEV V is given by. 

, ^ H 




(57) 
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where H is the magnetic held along the cr-direction. We also have 


Uijk = 4((5ij0fc + (cyclic)) f + 8(j)i(j)j(j)kf'' (59) 

Uijki = 4:(6ij6ki + (cyclic))f' + 8(5ij(j)k(l)i + (cyclic)) f' + (60) 

etc. Here /', /" etc. are the derivatives of / w.r.t. its argument. Since 




TTiTTj \ min 


= mlSii 


7r"0 


( 61 ) 


we get, 




2f(v^) 


H 

V 


(62) 


Also, 


ml = 2f(v^)+4v^r(v^) (63) 

It can be seen that Eqns. ^ and ^ together are equivalent to the k ^ 0 
limit of the identity in Eq.^. /^From Eq. 0 ,the IPI pion four - point vertex 
at zero momentum is given by 


f/4vr = 4:(6ij6ki + cyclic) f"(v^) (64) 


and from Eq.^, the crTTvr IPI vertex is given by 

= 4vS,,f(v^) 

Hence the tt — vr-scattering amplitude at zero momentum |]^ is, 

16u2(/"(x;2))2 


A47r = i(Sij6ki + cyclic..)[4:f" - 


2/'(u^) + 4u^/"(u^) 


(65) 


( 66 ) 


In other words 

+ (67) 

Irrespective of /", is at least as soft as H. This is equivalent to the soft 
pion theorems of current algebra 0. Note that this result has been proven 
here without the detailed assumptions of current algebra and PCAC. 
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On introducing the shorthand notation 

^ - ^in-lin CycHc] ( 68 ) 

and 

0*li2...in = (69) 

one can write the general result 

Ui^i 2 ..iN — ^ cyclic (70) 

k,l 

2k+l=N 

Now let us consider the zero momentum IPI vertex function Vnam-K (with 
N = n + m) corresponding to n-longitudinal and m-transverse external 
legs.Then all tensors will vanish when / > n as then at least some 

of the (j)'’s in these tensors will have to be transverse.Furthermore the leading 
IR behaviour is equivalent to treating v as being very large (compared to the 
momenta and equivalently to h ). It is clear that due to the exchange of the 
massless goldstone particles in the intermediate state, the n-point functions 
should have infrared singularities for sufficiently large n. From very gen¬ 
eral considerations, it follows that the nature of these singularities becomes 
more severe with increasing n. Thus, for sufficiently large no, the derivative 
j«.o(^ 2 ) structure 

= Eo(/Vt)-'-' + c„ + ^c,(/V"))' (71) 

fc>0 l>0 

where 

= - (72) 

V 

Therefore 

/(™+i)(ff2) = _ ^ + (73) 

A:>0 J 7 i>0 J j 

Now three situations can arise: (a) /^^cr^) is hnite or goes to zero slower 
than f'(cr^) as a ^ v. Then /("-o+i) is more singular than Also ^ 4 ,^ — 

H from (67): (b)/"((T^) and have the same leading infrared behaviour. 

Further, (63) already implies that ml and ^ which is a rather strong 

result. Then all for n > no have the same infrared behaviour. We shall 


20 



show shortly that this is inconsistent, (c) /”(<^^) vanishes faster than f'icr'^) 
as H —0. This has the consequence that has softer IR singularity than 
fin 2 ) fQj. ^ gy virtue of (63) this choice also implies that = 

This choice also implies that ^ 4 ^ vanishes faster than H from (67). As 
these features are already invalidated in perturbation theory, (c) will not be 
considered any more. 

Now, for both the situations a) and b) it follows that the dominant contri¬ 
bution to Tnarrnz will come from I = n {k = m/2) in (70); this is the maximum 
possible value for /. Hence 

r„ ~ (74) 

This further implies that 

^namn ^ ^{2n-\-m)TT (75) 


This means that as far as the infrared behaviour of all the vertex functions 
is concerned 

a ~ TT^ (76) 

This and the fact that pions are non-interacting in the IR leads to the stated 
result on xl viz. xl ^ H~ for 2 < d < 4 and xl ~ logH for d = 4. 

But these behaviours for xl are inconsistent with the situation (b) which 
demand xl ~ id“h Thus only (a) need be considered. 

In order to get an appreciation for the nature of the effective potentiaI,we 
analyse it in the large N limit.From the generating functional r[0. A] of Sec.3. 
one can recover the usual generating functional r[0] by eliminating A using 
its equation of motion. Since we are interested in constant helds 0, it is 
sufficient to have the A held also constant. Subtracting Eqn. ^ from the 
analog of Eqn. ^ that results from the renormalised Too (eqn |^), we get. 


(A-A„)|i- 


d^k _ 1 _ 

(27r)'^ (fc2 — A)(/c2 — Ao) 


(fc^ — /i)2 


)] 


Which can be solved to yield the implicit solution 

A = /(02-n2) 


(02-n^) (77) 

(78) 


for given Aq, m.T herefore the effective potential has the form. 


froo(0) 


4m 



(79) 
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where A is to be substituted as in Eqn. |^. Further we know that the min¬ 
imum of the effective potential is at, 0 ^ = v‘^,X = Aq such that Aq —>■ 0 ,u 7 ^ 
0, uAo = h, as h —>• O.This means that A and also Uoo{4>^) when expanded 
about the minimum 0^ = has Taylor coefficients which are increasingly 
singular in Aq and hence in h. 


6 Use of Schwinger - Dyson Equations 


In this section we present an argument based on Schwinger - Dyson equations 
that shows the diverging xl as a consequence of the infrared divergence of 
the Goldstone bosons. 

The heirarchy of Schwinger-Dyson equations are obtained by exploiting 
the following property of functional integrals 

0 = / (80) 
J d(j)i{x) 

To be specihc let us consider 

= j d^x 


5(0) 


[^ 0 ^ 0 - 0 ^ 0 -^/rV- 0 -^^A( 0 - 0 ) 2 ] (81) 


The above identity implies the following equation for T, the generating func¬ 
tional of IPI diagrams: 


6T 

64>i{x) 


where 


and 


+ D{x,k;x,k)(j)i{x) + 2D{x,k;x,i)(l)k{x) 
bi\ 

+D{x, k] y, j)V{x\ k'] y, j; z, l)D{x\ k'] x, k)D{z, /; x, i)] (82) 


D ^{x,a]y,h) 


50a(a;)(50b(|/) 


V{x,a-,y,b;z,c) 


6^T 

5(j)a{x)5(j)b{y)5(j)c{z) 


(83) 

(84) 
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It should be noted that D{x, a; y, h) and V{x, a; y, 6; z, c) are the propagator 
and 3-point vertex respectively only when the above derivatives are evaluated 
at the minimum of F namely, at a such that 

= 0 (85) 

d(pi 

Hence D and V in Eqs(|83[) and are 0-dependent. 

It will prove useful to display explicitly the behaviour of the IPI graphs in 
the large-N limit. For example, the four-point vertex has a 1/A^ dependence 
for large N. Likewise, the expectation value of the order parameter has a 
behaviour,and hence the 3-point functions in the ordered phase will have a 
factor. 

Now consider the SD equations for F^^. This equation is best repre¬ 
sented diagrammatically as in Fig. 2. In these hgures solid lines represent 



Fp2 

Figure 2: Schwinger-Dyson equations for the a two-point function 


exact (j-propagators and dashed lines exact 7r-propagators( except the lines 
with subscript ‘0’ denoting tree level propagators ), and all blobs represent 
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exact vertices. 


We assume some infrared regulator like the external field H. Then 
the Goldstone boson mass is not zero and we can safely study the zero mo¬ 
mentum behaviour of various graphs. The degree of infrared divergence is 
reflected by the leading m.n--dependence. Our input is that the exact pion 
propagator has the behaviour for k —>■ 0, i.e. we put in the infor¬ 

mation that the pion has a canonical behaviour in the infrared and study the 
consequences for the cr-propagator.Diagram 2(a) contributing to T^^ would 
behave like if To-tttt were 0(1) in the IR (and consequently the total T^^ 
would also be 0(1) due to the Ward identity (eqn ^)). This contribution is 
of order 0(1) in the large-N count. Thus, either the IR-singular behaviour 
of 2(a) is cancelled by other diagrams in £g (2) or To-TTTr must be at least 
as soft as Diagrams with potential IR-divergences are 2(b),2(c) and 

2(d); apart from all of them being of 0(1/N), none of them is as IR-singular 
as 2(a). Hence cancellation is ruled out and Tg-jr^ must behave at least like 
^4-d+Q a-positive semidehnite. Consequently T^^ must vanish as either 
or 777,4-'^+“ whichever is dominanat. 

While this proves conclusively that xl must diverge, we still need to restrict 
a further. It should be noted that T^^ has terms of order and 0(1) in 
the large-N count. On the other hand T^^ must vanish as some power of 
This is only possible if the terms are cancelled by such terms from 
the other diagrams. As all the diagrams except 2(a) are 0(1/A^), generically 
this can only happen from the ultraviolet contributions to 2(a). Within 1/N 
perturbation theory at least, such cancellations are not possible in which case 
the terms must be generated in the infrared by 2(a) and this is possible 
only if a = 0. This completes the proof based on SD-equations. 


7 Some Universal Features of Multipion Pro¬ 
cesses 

We have thus established that the divergence of the longitudinal susceptibility 
is very general depending only on the goldstone phenomena. In fact, we have 
shown on very general grounds that the leading infrared behaviour of the 
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longitudinal two-point function is completely universal. 

Now we show that the same considerations lead to certain universal fea¬ 
tures of multipion scattering amplitudes. Such universal features were al¬ 
luded to in Weinberg’s article 0. Consider, for example, the expression (67) 
for the vr —TT scattering amplitude at zero momentum. Eqn (74) can be recast 
as 


/' ^ 2 < d < 4 

~ {lnH)~^ for d = A 

Hence 4u^/'u^) » H as id —>■ 0 and one has 


( 86 ) 

(87) 



( 88 ) 


On recalling that in the absence of the infrared cut-off H, the scale that 
replaces H is E^, we recognise that the leading universal infrared behaviour 
^ 4 ^ ~ E‘^ is indeed the well known result for soft pion scattering (usually 
derived in d = 4). This universal behaviour does not depend on the space 
time dimensionality. 

However (^) demonstrates that even the next to leading term is universal. 
The above considerations can be applied to multi-pion scatering amplitudes 
and one again Ends the same universality of the hrst two terms. It is also 
clear from this paper that these results are non-perturbative in nature. 

The presence of these non-analytic (in energy-momentum variables) terms 
in d = 4 had been noted by Li and Pagels 0, the coefficients of which in 
the case of tt — tt scattering had been calculated by Lehmann, and Lehmann 
and Trute 0 using unitarity. Weinberg showed that the application of renor- 
malisation group can be used to show the universality of these logarithmic 
terms. He implicitly uses chiral symmetry. Though superhcially it appears 
that the universality of the logarithmic terms is a one-loop result, careful re¬ 
examination of Weinberg’s demonstartion shows that it is an all-loop result. 

In our approach, the origin of the universality of these subdominant log¬ 
arithmic terms is directly related to universal divergence of the longitudinal 
susceptibility. 

Our large-N derivation, which links the divergence of the xl to the special 
properties of the auxiliary held also demonstrates the universal features in a 
very straightforward manner. 
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The multipion scattering amplitudes in this case can be classified into 
three categories: a) those involving only A and vr propagators, and the Avr^ 
vertices. For the scattering amplitude involving 2n pions (including initial 
and final states ),one needs n — 1 A and n — 2 tt propagators. Hence the 
IR-behaviour of all multipion amplitudes of this class is ~ b) Amplitudes 
involving one nX vertex and n Att^ vertices. The nA-vertex scales like 
Hence this contribution to the amplitude scales like (p^)'^'^^, which is less 
dominant than the behaviour a). Finally, there is the class c) which consists 
of amplitudes involving combinations of multi-A vertices with m < n.Their 
contribution can be seen to be even less dominant than b). 

Because of the k^I{k^) factors in the denominators of the A-propagator, 
one gets a universal correction factor and the leading IR behaviour 

of these correction factors is the same in the linear and non-linear a-models. 
This is precisely the same quantity that governs the universal divergence of 
XL- 

We now make a few geenral remarks about conventional chiral perturba¬ 
tion theory (ChPT) and compare it with our approach. 

As we have already stated many times earlier, the leading IR behaviour 
of systems with Goldstone bosons is completely universal, and in the sector 
of Goldstone bosons only, is represented by the Lagrangian 


C 




1 (tt X 

J TT 


This, following Weinberg, can be cast in the form 


(89) 


( 90 ) 


The validity of (l89|) is, however, strictly in the extreme infrared. One of 
the procedures adopted to extend to regions beyond the extreme IR is 
the so-called chiral perturbation theory. In the form that is mostly practised, 
GhPT amounts to starting, for example, from a Lagrangian of the form 


C = Cr 






( 2 ) 


{D,'K-D^TyY + ^-X-{D,Ty-D,'K){D^Ty-D,'K) + - ■ ■ (91) 


Gontrary to the attitude of the early days of effective Lagrangians, one 
now carries out loop calculations with Lagrangians of the type (|90|) and (pl|). 
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Owing to perturbatie non-renormalisability of these lagrangians, one needs 
infinitely many counter terms to absorb the various resulting inhnities.. How¬ 
ever, to any any desired accuracy in the momenta involved, one needs only a 
hnite number of counter terms. The structure of the resulting renormalised 
Lagrangian whose domain of validity is somewhat larger that that of is 
that it is charaterised by a number of unknown parameters corresponding to 
the number of counteer terms used as well as certain universal parts inde¬ 
pendent of these unknown coefficients (/jr is always treated as known). This 
is indeed the general structure argued by Weinberg. 

Our analysis that the leading IR behaviour is governed by the Gaussian 
fixed point Qn-i points to a more efficient way of realising the results of 
ChPT. This suggestion was implicit in also. This consists of perturbing 
the IR fixed point directly with the help of irrelevant operators. Their influ¬ 
ence on all correlation functions can be directly computed by the use of the 
RG techniques. The ultraviolet degrees of freedom have no role to play in 
this scheme. This method will be applied to practical calculations of GhPT 
elsewhere. 


8 Conclusions 

In this paper we have emphasised that the quantum fluctuations of the 
Goldstone bosons have some drastic effects which are not encountered ei¬ 
ther in a tree level or loopwise analysis of the linear or the non-linear a- 
models.Nevertheless,the leading infrared behaviour of the correlation func¬ 
tions again become identical in the linear and the non-linear versions.lt is 
remarkable that the longitudinal fluctuations also become ‘massless’ as a 
consequence of the masslessness of the Goldstone bosons. We have exhibited 
the generality of this phenomenon giving three separate arguments based on 
i). Renormalization group flows, ii). Ward identities, iii). Schwinger-Dyson 
equations. 

Heuristic derivations of the soft - pion theorems within the linear cx-model 
are based on the fact that a-mass (xl^) being non-zero, for all processes 
involving momenta much smaller than this mass the cx-propagator can be 
replaced by One may wonder whether the results of this paper would 

vitiate this and consequently the soft-pion theorems too. But as demon¬ 
strated in the sections on the large-N analysis as well as the one on Ward 
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identities, the soft pion theorems are still valid. The deeper reason for this 
is, of course, the renormalisation group point of view according to which the 
leading IR-behaviour is controlled by the Gaussian hxed point. 

As stressed sufficiently in this article, both the qualitative and quanti¬ 
tative aspects of the result,as far as the susceptibilities are concerned, were 
known in the literature from somewhat differing perspectives which were not 
general enough. What we have done here is to demonstrate the full generality 
of this result as well as extend the analysis to include the nonperturbative 
IR behaviour of correlation functions also. Despite the fact that the result 
concerning xl had appeared in the literature in different guises , its signih- 
cance does not appear to have been well appreciated. 

We emphasise that in the chiral limit the ‘sigma particle’ has a singular 
infrared behaviour (i.e. logarithmic divergence) even in four dimensions. 

Independently of the space-time dimensions, the Goldstone bosons have 
canonical dimensions and decouple from each other and their logitudinal 
partners in the infrared.We pointed out that this can be understood from the 
Goldstone theorem, current algebra and PGAG on the one hand and the flow 
towards 0{N — 1) Gaussian hxed point on the other hand. This latter point 
of view is very valuable because one can easily study the non-leading correc¬ 
tions to the infrared behaviour by standard renormalization group methods 
[§]. All effects of the microscopic theory could be parametrised in a few irrel¬ 
evant and marginal operators around the 0{N — 1) Gaussian hxed point and 
their ehects computed. This approach provides a technically more efficient 
and conceptually simpler alternative to chiral perturbation theory, which 
focusses on the infrared divergence of the Goldstone bosons instead of the 
non-renormalizable ultraviolet divergences. 

We showed that the leading N approximation takes into account all the 
infrared behaviours accuratly, including that of the longitudinal huctuations, 
in contrast to perturbation theory based on the linear or the non-linear a- 
models.We developed a technique which allows explicit computation of the 
correlation functions even far away from the infrared region.This provides a 
useful non-perturbative technique for obtaining a qualitative understanding 
of the effects of ffuctuations. 

A convenient tool for understanding longitudinal susceptibility is the 
equation of state p.In terms of the reduced variables y = H/M^^x = 


where M is the magnetisation, t the deviation from the critical tem¬ 
perature,and 6, /3 the standard exponents, the equation of state takes the 
universal form 


y = fix) (92) 

As H tends to zero, x must tend to one of the roots of fix). With suitable 
normalisation this root (corresponding to spontaneous magnetization) can 
be arranged to be at x = —l.Now,the question of whether xl diverges or not 
can be answered by knowing the manner in which f{x) approaches zero as x 
approaches — l.If this approach is parametrised as 

fix) = il + x)^^P (93) 

where p is a positive semidehnite number, the behaviour of xl as H tends to 
zero is given by 


XL ~ (94) 

In the large-N expansion for arbitrary d (except d=4 where there are loga- 
rthmic modi£cations),the equation of state is 


y=il + x) 


2 

d-2 


(95) 


as can be proven using the euclidean version of eqn(^) of sec 4(see also 
[|]). As the exponent of (1 -|- x) is greater than 1 one concludes that the 
longitudinal susceptibility diverges for all T < Tc. Modulo the wavefunction 
renormalisation constant, Xl is just 

There are interesting consequences of this phenomena both for finite tem¬ 
perature QCD as well as QCD at zero temperature. In the case of the former, 
under the well argued scenario that this is a second order transition, Wilczek 
and Rajagopal conclude that the transition must belong to the universal¬ 
ity class of d = 3 0(4) magnet models. Then the general results of this 

paper become applicable. 

In numerical simulations of the QCD chiral phase transition at hnite tem- 
perature,the quantity that has been used to characterise the nature of the 
phase transition is the so-called A-cumulant dehned by 


A = 


dlmp'ip 

dlnrrin 


( 96 ) 
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where V’V’ is the chiral condensate and niq the current-quark mass.Thus 

T > Tc 'ijj'ijj ^ mq A = 1 

T = T, ~ A = 1/^ (97) 

Below Tc,but close to it, one should expect, on the basis of the considerations 
of this paper 

T <Tc A ~ mjA (98) 

The same considerations apply to the spontaneously broken phase of zero 
temperature QCD. Here one should hnd 


A ~ mqlnrriq 


(99) 


characterstic of Goldstone phenomena in four dimensions. 

Of course, the arguments of Wilczek and Rajagopal that maps the problem 
to a classical stat mech problem in d = 3 is valid exactly at T^. But continu¬ 
ity would demand that at temperatures in the vicinity of Tc too the effective 
dimensionality should be close to 3. Thus, the expectations based on eqs 
and ^ ) are the correct ones. This should yield the experimentally 


interesting signal that for QCD at hnite temperatures T < Tc, xl should 


smoothly go over from eqn (pSD to eqn(^) as T is lowered. 

Indeed, wherever there are Goldstone bosons, one should expect to see the 
singular behaviour of xl- It is therefore somewhat disappointing that clear 
experimental signatures of this effect have not yet been found even for the 
classic ferromagnets.The lack of isotropy in most real life ferromagnets would 
make the experimental establishment of this effect a challenging one. An at¬ 
tempt to see this for ferromagnets has recently been made by [^. They 
should have been looking for a xl ~ behaviour rather than xl ~ 

^-1/3 rjpg behaviour of xl is ’super-universal’ in the sense 

that it depends only on dimension and the existence of Goldstone modes 
should be established in as many diverse experimental circumstances as pos¬ 
sible. As Patashinsky and Pokrovsky |T^ have pointed out, the effect should 
also be observable in the superfluid phase of He.Numerical simulations of the 
0{N)- models in three dimensions, performed in the ordered phase, should 
also establish this effect. 

Kogut et al. report a completely different behaviour and different expo¬ 
nents for the hnite temperature chiral transition in QCD in comparison to 
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[Q], In particular, Kogut et al. do not find the crucial softening of the longi¬ 
tudinal mode. More speci£cally,they argue that PCAC restricts the equation 
of state to be of the form 


f{x) = {l + x) (100) 

There is a fallacy in their argument; PCAC implies the above form of the 
equation of state if and only if the longitudinal susceptibility is assumed to 
be hnite. Only then a Taylor expansion around zero applied magnetic held 
is permissible.But our results about the divergence of xl imply that the as¬ 
sumption of a hnite xl is incorrect. Once allowance is made for a diverging 
Xli PCAC no longer places any restriction on the form of the equation of 
state. Also, Kogut et al. carry out their analysis in d = 4 even at the critical 
temperature. 
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